Topological pumping in the one-dimensional Bose-Hubbard model 



o 

(N 
D 

(N 



I 



I 

o 
o 



> 
m 
rn 

O 

cn 



Davidc Rossini, Marco Gibcrtini, Vittorio Giovannetti, and Rosario Fazio 
NEST, Scuola Normale Superiore and IsUtuto Nanoscienze-CNR, 1-56126 Pisa, Italy 

By means of time-dependent density matrix renormalization group calculations we study topo- 
logical quantum pumping in a strongly interacting system. The system under consideration is 
described by the Hamiltonian of a one-dimensional extended Bose-Hubbard model in the presence 
of a correlated hopping which breaks lattice inversion symmetry. This model has been predicted to 
support topological pumping [E. Berg, M. Levin, and E. Altman, Phys. Rev. Lett. 106, 110405 
(2011)]. The pumped charge is quantized and of topological nature. We provide a detailed analysis 
of the finite-size-scaling behavior of the pumped charge and its deviations from the quantized value. 
Furthermore we also analyze the non-adiabatic corrections due to the finite frequency of the modu- 
lation. We consider two configurations: a closed ring where the time-dependence of the parameter 
induces a circulating current, and a finite open-ended chain where particles are dragged from one 
edge to the opposite edge, due to the pumping mechanism induced by the bulk. 

PACS numbers: 05.30.Jp, 64.70.Tg, 67.85. De 



I. INTRODUCTION 

A charge/spin flow can take place, in the absence of 
any direct bias, by a periodical modulation of some pa- 
rameters (magnetic fields, coupling constants, external 
gates, . . . ) affecting the dynamics of the system^. This 
transport mechanism, known as quantum pumping, is 
intimately related to the properties of the instantaneous 
eigenvalues and eigenfunctions of the system Hamilto- 
nian. Pumping is adiabatic if the variation of the param- 
eters is slow as compared to the characteristic time scales 
of the system. Since the pioneering work of Thouless in 
1983^ pumping has been attracting an intense theoretical 
interest^''^ and experimental interesf*^. In the scattering 
approach to quantum transport^, the pumped charge in 
an adiabatic cycle can be expressed in terms of an integral 
in parameter space which does not depend on the evolu- 
tion period T. Accordingly, adiabatic quantum pumping 
appears to possess an intrinsic geometric character which 
allows us to relate it with the Berry phase accumulated 
by the system wave function during the cyclic evolution^. 
Interestingly enough, under proper conditions, the quan- 
tum pumping can acquire also a topological nature which 
yields to a quantization of the pumped charge/spin per 
cycle. This effect is connected to the appearance, in the 
Hamiltonian spectrum, of gapless points enclosed by the 
adiabatic loop in parameter space^'^. Conditions for the 
quantization of the pumped charge have been discussed 
in Ref. 7, while, for the case of spin transport, a clas- 
sification of topological non-interacting pumps has been 
given in Ref. 8. 

The study of quantum pumping has been extended 
also to interacting systems^"^^. Here the Hamiltonian 
gapless points responsible for topological pumping stem 
from many-body effects and are typically associated to 
quantum phase transitions. This suggests that quantized 
adiabatic transport can be achieved by evolving the sys- 
tem around a loop in parameter space enclosing a criti- 
cal point. Recently it has been suggested, with the sup- 
port of an analysis based on bosonization, that quantized 



pumping occurs in systems of interacting lattice bosons 
at integer filling in one dimension^'^. When both on-site 
and nearest-neighbor interactions are taken into account, 
a critical point arises separating two distinct insulating 
phases^^. Encircling this point in parameter space adi- 
abatically is expected to pump exactly one boson per 
cycle. 

The interesting result given in Ref. 13 calls for a de- 
tailed analysis of the pumping mechanism in the ex- 
tended Bose-Hubbard model. Here we tackle this prob- 
lem by means of numerical time-dependent density ma- 
trix renormalization group (t-DMRG) simulations^'^. We 
compute the number of pumped particles per cycle as 
a function of the system size and show that the scal- 
ing is consistent with unit quantized pumping when the 
adiabatic loop encloses the critical point, while it oth- 
erwise vanishes. We show that the evolution period T 
does not affect the pumping properties of the system, 
provided it is chosen so that the system remains adia- 
batic with respect to the minimum instantaneous gap. 
The topological character of the pumping mechanism is 
tested by considering different loops and assessing the 
independence of the number of pumped particles on the 
specific trajectory in parameter space. A quantitative 
analysis of finite-size and non-adiabatic effects is a neces- 
sary and important step towards a possible experimental 
observation of topological pumps. 

Pumping in the Bose-Hubbard model can be ana- 
lyzed in two different configurations. One can consider 
a ring geometry. Here the modulation of parameters 
will result in a persistent pumped current. Alterna- 
tively one can place the interacting region between two 
(non-interacting) reservoirs, in which case particles are 
pumped from one lead to the other. We will consider in 
great details the first case and comment on the second 
setup at the end of the paper. Concerning the latter case, 
we will provide results about the simulation of a pure 
interacting system in an open-ended geometry, showing 
that the ultimate effect of the pumping mechanism after 
one loop is to drag particles from one edge to the opposite 
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edge of the chain, leaving the bulk unaffected. 

The paper is organized as follows. In the next section 
we define the model and summarize the salient features of 
its phase diagram which are of relevance for the topolog- 
ical pumping. Section III is devoted to the description of 
the method and of all the important technical details we 
employed to study both static and dynamical properties 
of the model. We then discuss quantitatively the behav- 
ior of the gaps at the critical point and along the loops 
in parameter space we are going to consider (Sec. IV). In 
Sec. V we explain how we compute the pumped charge. 
We first consider the ring setup (Sec. VI), discussing the 
deviations from perfect quantization due to the finite 
number of sites and the finite modulation period. Then 
we address the case of a transport experiment in an open- 
ended interacting system, where particles arc pumped 
from one end of the chain to the other (Sec. VII). The 
last section is devoted to the conclusions of this work. 



II. THE MODEL 

The model we consider consists of a collection of inter- 
acting lattice bosons on a chain of L sites at unit filling 
(i.e., with the total number of particles set to be con- 
stant and equal to the number of sites). For this system 
we assume an Hamiltonian given by the sum of two con- 
tributions: 
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FIG. 1: Phase diagram of the one-dimensional extended Bose- 
Hubbard model defined by the Hamiltonian ?^o, at unit fill- 
ing^^. The various regions separated by lines denote the fol- 
lowing phases: MI = Mott insulator, HI = Haldane insula- 
tor, SF = superfluid, DW = density wave. The boundaries 
between the three insulating phases have been obtained by 
computing the string Cstring and the staggering order param- 
eter Odw, while the SF phase has been located by analyzing 
the closure of the ground-state charge gap (see Ref. 18 for 
details). In the whole paper we concentrate on the region at 
[7 = 5 and close to the critical point C at Vc 2.7 (filled 
circle), as indicated by an horizontal dotted line. The values 
of U and V are expressed in units of the hopping strength. 



H = no + nx . 



(1) 



The first term describes an extended Bose-Hubbard 
model, defined by the expression 
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where, for j — 1, ■ ■ ■ , L, the operator bj (bj) creates (de- 
stroys) a boson at lattice site j [hj = Sjfej being the as- 
sociated occupation operator), while, depending on the 
chosen boundary conditions, ^l+i is set equal to hi (ring 
geometry) or to (open chain geometry). In Eq. (2) 
U and V respectively denote the on-site and nearest- 
neighbor interaction strengths in units of the hopping 
strength, which hereafter implicitly sets the energy scale. 

The zero-temperature phase diagram of 'Ho in the 
U — V plane is sketched in Fig. 1. For small values of 
the nearest-neighbor repulsion there is a direct super- 
fluid (SF) to Mott-insulator (MI) transition. For suf- 
ficiently large values of J7, before entering a density- 
wave (DW) phase characterized by a non-uniform boson 
density with staggered ordering^^, the model defined by 
Eq. (2) exhibits a second order quantum phase transi- 
tion between two insulating phases^^'^^'^®, the distinc- 
tion among them being protected by the bond-centered 
inversion symmetry of 'Ho^^'^^''^^ . Specifically, upon in- 
creasing V the system evolves from a conventional MI 
to a bosonic Haldane insulator (HI), which displays 



similar correlations to the Haldane gapped phase^^ of 
spin-1 chains (it is indeed identified by a non- vanishing 
string order parameter^'', defined according to ©string = 

where 5n^ denotes the 
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boson number fluctuations from the average filling) . No- 
tice however that while for the ring configuration, i.e., 
for Periodic Boundary Conditions (PBC), the HI phase 
is associated to a single non degenerate ground state, 
for Open Boundary Conditions (OBC) it exhibits a four- 
fold dcgeneracy^^ that needs to be accounted for, when 
designing the parameter loop (sec Sec. VII for details). 
The MI-HI transition defines the region of interest for our 
analysis: in the rest of the paper we will fix the on-site 
interaction to the value ?7 = 5, so that the associated 
critical point C occurs at Vc ~ 2.7 (see filled circle and 
dotted line in Fig. 1). 

As discussed in Ref. 17, in order to have a non-zero 
pumped charge, one should drag the system adiabati- 
cally along a loop F enclosing C in the parameter space 
of the Hamiltonian. For this purpose, a second contribu- 
tion T-Lx is added to Eq. (2), which allows for an adiabatic 
connection between the MI and HI phases, by explicitly 
breaking the lattice inversion symmetry of "Hq. In partic- 
ular, wc choose a correlated hopping between neighboring 
sites which reads 
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FIG. 2: (color online). Parameter space of the Hamiltonian 
in Eq. (1) at (7 = 5 and unit filling. The point (Vc,0) with 
K; = 2.7 is an isolated critical point separating two half lines, 
iy < Vc,0) and {V > Vc,0), along which the system is in 
the MI and in the HI phase, respectively [the same clear-cut 
distinction between the two phases being not possible when 
A 7^ 0, owing to a lack of inversion symmetry of Eq. (3)]. 
Different colors denote the loops F that we have considered 
for the adiabatic evolution of the system. The loops marked 
with black circles {"Loop 1"), red squares ["Loop 2"), and 
green diamonds ( "Loop 3") enclose the critical point, and are 
thus expected to be associated with a quantized number of 
pumped particles. The loop with blue triangles ( "Loop 4 ") is 
instead expected to give null pumped charge, as it does not 
encircle (Vc, 0). 

with A being a (controllable) intensity parameter. With 
this choice, upon fixing U as detailed above and by exter- 
nally controlling V and A, we can now drag the system 
along a pumping cycle F around the isolated critical point 
(y, A) — {Vc, 0) that represents C in the X — V plane (see 
Fig. 2). Specifically, assuming the system to be initial- 
ized into the ground state {ipcs) of the Hamiltonian (1) 
evaluated for some given point (Vo, Aq) of F, we evolve it 
via a time-dependent Hamiltonian 'H[V{t), X{t)] obtained 
from (1) by varying parametrically in time the quantities 
V and A along the loop F (the trajectory being periodic 
with period T). Accordingly, at time t the vector state 
becomes 

|^(0)=^e-^o«[nO.A(0]d*'|^^^^^ (4) 

with {V{t),X{t)) e F for all t and with (F(0),A(0)) = 
{yo, Xo)- This is the state on which, as detailed in Sec. V, 
we are going to evaluate the pumped charge accumulated 
on a cycle. 

III. METHODS 

In order to study the many-body system in Eq. (1), 
we employ a numerical approach based on the t-DMRG 
method in the formalism of the matrix product states 
(MPS)^^. In particular, to characterize the properties 
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FIG. 3: (color online). Implementation of PBC using a OBC 
scheme for the underlying MPS ansatz. The left picture 
exemplifies the standard approach, in which each physical 
site (black circles) is associated to a computational site (red 
boxes), so that the MPS has a cyclic structure. The right 
picture denotes our scheme: to each physical site correspond 
two computational sites, thus effectively "flattening" the ring 
as an open-ended ID shape. In the drawings we sketched a 
chain with L = 8 sites. 



of the ground state iV'Gs) of Eq. (1) (such as the energy 
and the relevant correlation functions), we impose a MPS 
trial wave function, and perform a variational minimiza- 
tion of the energy cost function site by site. Subsequently, 
in order to follow the evolution of the vector in Eq. (4), we 
adopt the time-evolving block decimation (TEBD) pro- 
cedure^'^. In our simulations we were able to keep track 
of the globally conserved number of bosons, thus always 
working in the canonical ensemble. We fixed the num- 
ber of particles in order to have unit filling {h = 1), and 
distributed them over lattices of up to L = 200 sites^^. 

The choice of the boundary conditions of the chain re- 
quires a special attention: as a matter of fact, most of 
the time-dependent simulations presented in this paper 
(except those in Sec. VII) have been carried out for a ring 
geometry, i.e., under PBC. It turns out that the standard 
t-DMRG algorithm is not optimized for these boundary 
configurations. Indeed, on one hand this choice imposes 
a cyclic structure to the MPS ansatz^^, thus resulting in 
a considerable slowdown of the static optimization algo- 
rithm, unless some clever strategies are adopted^^. On 
the other hand, the inability to construct an orthonormal 
basis, starting from the leftmost and from the rightmost 
MPS sites and requiring the proper unitary conditions 
on the defining matrices, prevents from the possibility 
to write any bipartition of the system in two contigu- 
ous blocks of sites as a Schmidt decomposition. This 
causes a failure in the standard TEBD scheme and leads 
to an uncontrollable truncation error, which demands for 
a substantial increase of the required MPS bond-link di- 
mension m and a drastic slowdown of the computational 
performances. 

Instead of directly imposing PBC, wc found it more 
convenient to implement a standard t-DMRG scheme 
with OBC, where the ID physical system with a ring 
geometry has been "flattened" as a narrow string and 
therefore mapped into an open ended chain. As pictori- 
ally shown in Fig. 3, we associate to each computational 
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site two distant physical sites. This trick enables us to use 
a conventional t-DMRG code for open boundaries and to 
take advantage of all the above mentioned methods for 
standard OBC. In particular we are able to easily imple- 
ment the conservation of abelian symmetries and to con- 
trol the truncation error within the TEBD scheme. The 
price we have to pay is that of doubling the local Hilbert 
space of each computational site, and using a larger m 
value in order to get accuracies comparable to the OBC 
situation^'^. 

Our code dynamically updates the bond-link dimen- 
sion used at each link, both in the static algorithm and 
during a single step of time evolution, according to a 
given threshold value in the sum of the discarded weights 
£. We found that, in order to keep it below 5 x 10~^, a 
number of states m < 400 is generally sufficients^. We 
checked that a further reduction of e does not signifi- 
cantly affect our results. In the TEBD dynamics, we 
employed a fourth-order Trotter expansion of the time- 
evolution operator, thus minimizing the error due to a 
finite time step (but still admitting time steps of the or- 
der (it - 5 X 10-2). 



IV. ENERGY GAPS 

As we are dealing with finite-size systems character- 
ized by the Hamiltonian in Eq. (1), the behavior of the 
ground-state energy gap is important, especially at the 
critical point and along the loop in parameter space. In 
order to enter the adiabatic regime, ideally one would like 
to complete the loop in a time T which is much larger 
than the inverse of 



As„i = min Ay,A 
(KA)er 



(5) 



that is the smallest of the gaps Ay. a of the Hamiltonians 
H[y, A] evaluated along the selected loop F. In addition, 
for a finite-size system, the gap Ac of the Hamiltonian (1) 
at the critical point [U = 5, 14 — 2.7, A = 0) is never ex- 
actly zero. Therefore, on decreasing the amplitude of 
the loop enclosing the singularity, one should start to ob- 
serve deviations from the quantized value for the smallest 
loops for which the finiteness of the critical gap starts to 
be resolved. An additional constraint on the amplitude 
of the pumping path raises due to the emergence of the 
DW phase close to the Haldane insulator. 

All these requirements lead us to consider the loops 
shown in Fig. 2, where we depict four paths in the 
\ — V parameter space along which we have computed the 
evolution governed by the time-dependent Hamiltonian 
-Hit) = -H[V{t),\{t)] with period T: Hit + T) = n{t). 
Three of them (black, red, and green set) encircle the 
critical point (14,0), while the last one (blue set) is triv- 
ial, meaning that the gap is never vanishing inside it. 
The asymmetric shape of the loops with respect to the 
critical point is a trade off between the need to avoid the 
influence of the nearby DW phase, and at the same time 




FIG. 4: (color online). Finite-size scaling of the charge gap 
Av,\ at A = 0, close to the critical point and for U = 5. 
At the phase transition, 14 ~ 2.7, the gap scales to zero as 
L^^ (see the inset, where we plot the minimum reached value 
Ac as a function of L — the dashed line indicates a behavior 
Ac ~ 1/L and is plotted as a guide to the eye). Away from 
criticality, the gap is almost independent on the length, as 
long as L > 100. As detailed in text, the values of Av,\=o 
and Ac reported here have been computed by evaluating the 
charge gap Achargo of the system. 



to have a contour sufficiently far from the critical point 
in order to minimize the effects of the finiteness of the 
critical gap Ac. 

In order to define the time scales on a more quantita- 
tive basis, in the following we are going to analyze the 
scaling of the gaps Ac and A^a as a function of the sys- 
tem size L. To do so, we focus on the ground-state charge 
gap, that is the difference in the energy needed to add 
(A£'+) or to remove {AE~) a single boson in the system: 
Achargc = AE~^ — AE~ . In the whole MI, as well as in the 
HI phase far from the transition with the DW phase^"', 
the lowest energy excitations of the Hamiltonian (1) are 
indeed charge excitations, therefore the values of Achargo 
(evaluated for the proper V and A values) can be used to 
derive Ac and A^a- 

Let us first discuss the gap close to the critical point, 
Ac. In Fig. 4 we plot the size dependence of the charge 
gap, and show that, in the thermodynamic limit, it van- 
ishes as 1/L. For our convenience we computed the gaps 
in a chain with OBC. The differences are not relevant for 
our purposes (i.e., to determine the relevant time-scales 
for the pumping period). On the other hand in this way 
we could easily analyze the critical gap for chains up to 
few hundred sites^^. 

Now the interesting point which permits us to address 
the regime of topological pumping is that, already for 
L ^ 200, the gap Ac turns out to be almost one order of 
magnitude smaller than the typical gap A^.a encountered 
along the loops F of Fig. 2. A quantitative analysis of 
the value of the gap along one particular loop is shown in 
Fig. 5, where we took data directly in the case with PBC 
(specifically, here we analyzed "Loop 1"). In accordance 
with the fact that along the loop the system remains 
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FIG. 5: (color online). Finite-size scaling of the gap /S.v,\ 
along "Loop 1" (see Fig. 2) for three different system sizes, 
L = 60, 100, 140 [as in Fig. 4, the values of Av,a have been ob- 
tained by computing the charge gap Achargo]. The minimum 
value of the gap Asm (indicated by an arrow in the figure) is 
the bottleneck to define the adiabatic dynamics. 



FIG. 6: (color online). Excess energy, with respect to the 
actual adiabatic ground state, evaluated at the end of the 
loop, i.e., once the system parameters have been restored to 
the initial values. Here T denotes the total evolution time. 
The data have been obtained by performing a cyclic variation 
of the parameters (V^, A) along "Loop 1", starting from the 
upper point in the A — V plane. 



distant from criticality, the gap barely depends on the 
size. We notice that the value of Ay.x is not uniform 
along the loop and can significantly diminish at some 
points. This imposes strict conditions on the evolution 
period T necessary to observe quantized pumping and 
calls for a careful analysis of the role of non-adiabatic 
corrections. 



A. Testing the adiabatic condition 

Roughly speaking, one could expect that the adiabatic 
regime is achieved for periods T ^^m- ^ good es- 
timator of the quality of the adiabatic approximation 
along the performed loop is given by the excess energy 
Ai?oxc(0- This is quantified by the difference between 
the energy of the wavefunction during the evolution, 
{^{t)\nit)\iP{t)), and the energy (^gs (OI^WIV-GS W) 
corresponding to the instantaneous ground state |?/'gs(0) 
of 'H{t), i.e., 

Ai^cxcW = m)\nt)\m) - (^GsWi?iwi^GsW) ■ 

(6) 

In Fig. 6 we analyzed the value AEe^dT) at the end of 
the loop, as a function of the total evolution time T. The 
excess energy is generally decreasing with T, thus mean- 
ing that the adiabaticity condition along the loop can be 
progressively approached by slowing down the variation 
of the system parameters {X,V). Moreover, we note that 
Ai?exc(r) is not strongly affected by the system size L. 
This is consistent with the fact that the instantaneous 
gap barely depends on L, provided this is sufficiently 
large (see Fig. 5). We finally point out that, far from 
the adiabaticity condition, such energy is not guaranteed 
to be monotonic in T. 



V. PUMPED CHARGE 

Let us now introduce the quantity of interest, the 
pumped charge, to be evaluated by means of direct 
TEBD numerical approach for the time evolution of ID 
many-body systems. 

The number of particles which have flown through site 
j between time t and t' > t can be expressed as the 
integral of the instantaneous local bosonic current J7j (t) 
at site j: 

Q'f}. = J,{r) dr , (7) 

with Jj{t) ~ {il]{t)\Jj{t)\ip{t)) being evaluated as the 
expectation value of the local current operator Jj{t). 
Notice that, even though Eq. (7) is expressed in the 
Schrodinger picture, the operator Jj{t) depends on t via 
the parametric temporal dependence of V and A enter- 
ing the system Hamiltonian. Explicitly, such operator is 
determined via the continuity equation 

j,{t)-j,^^{t)^i[h„n{t)\, (8) 

yielding an expression, in terms of two-point correlators 
on nearest neighbor sites, of the form 

Jj{t) = ^( - h]bj+i + Xhfb]bj+i - h.c.) . (9) 

When integrating over a period T, the quantity (7) 
gives the number of pumped bosons per cycle. Under 
the conditions of topological pumping, this quantity is 
independent of the initial time t and of the starting point 
inside the loop, i.e., 

Q^^It = Qo,T ^ Qj (T) = i Jj (r) dT , (10) 

J loop 
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furthermore, owing to the equivalence between lattice 
sites, for a system with PBC, Qj{T) does not depend 
on j. 

In evaluating the function (10), we are going to con- 
sider two different configurations: a ring with periodic 
boundary conditions, and a finite chain with open bound- 
ary conditions. In the first one (Sec. VI), the periodic 
modulation of the external parameters {V{t),\{t)) in- 
duces a uniform and persistent current J{t) that is in- 
dependent of the position j. The pumped charge Q is 
the current passing through a link integrated over a pe- 
riod T. In the second setup (Sec. VII), the boundaries 
break the translational invariance and the instantaneous 
current Jj{t) is actually a site-dependent quantity. As 
we will show below, it turns out that the pumping mech- 
anism tends to deplete the particle density on one side, 
and to increase it on the opposite one. This progressively 
destroys the insulator, which can be stabilized only at a 
commensurate filling, starting from the edges. We mea- 
sure the current, and hence the pumped charge, close to 
the middle part of the chain, where the filling remains 
close to unity for the whole duration of the cycle. 



VI. RING CONFIGURATION 

As already mentioned in the previous section, for a sys- 
tem with PBC, the pumped charge on a cycle, Eq. (10), 
does not depend on j, i.e., Qj{T) = Q{T). As a matter 
of fact, the same property applies for all the quantities 
defined in Eq. (7), yielding Q'fl = Qt,t'- 



A. Finite-size scaling of the pumped charge 

Our numerical results for Q{T) after one cycle with 
a sufficiently long period T, in the sense discussed in 
the following Subscc. VIE, and enclosing an isolated 
quantum critical point are summarized in Fig. 7. There 
we report data (depicted by black circles, red squares 
and green diamonds) corresponding to the three differ- 
ent loops around the critical point in Fig. 2, as a function 
of the system size L^^. We note that the total pumped 
charge approaches the unit value for sufficiently large sys- 
tem sizes, no matter how we choose the specific path 
encircling the critical point. This is an important obser- 
vation supporting the topological origin of the quantiza- 
tion of Q, which, strictly speaking, is valid only in the 
thermodynamic limit. Remarkably, it rules out the possi- 
bility that Q remains simply a geometric quantity, even 
in the limit of an infinite chain, while instead it is has 
to be taken as a topological and quantized property of 
the system. As long as the size L of the system is finite, 
we have discrepancies in the final value Q at the end of 
the cycle computed for different paths. Such discrepancy 
decreases with increasing number of lattice sites and van- 
ishes in the thermodynamic limit, when Q is quantized 
to unity, independently of the specific path considered. 
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FIG. 7; (color online). Ring configuration: number of 
pumped particles Q per cycle as a function of system size 
L, corresponding to the three different loops encircling the 
critical point (black circles stand for "Loop 1", red squares 
for "Loop 2", while green diamonds for "Loop 3" - see 
Fig. 2). In our simulations we used different periods T ac- 
cording to the size, such to ensure the adiabatic approxi- 
mation along the whole loop {e.g., we took T — 100 for 
L = 8, and T = 350 for L = 140). Data are displayed in 
a logarithmic plot. The straight lines are best fits to the 
numerical data for L > 60, and are respectively given by 
1 - Q ?a 1.64 X e" ^ 
{"Loop 2"), 1-Q 



-0.031 L ( uj^^^p j„s^^ 1 - Q 1.76 X e" 
1.32 X e-^o^si {"Loop 3"). 



A deeper investigation of the scaling properties of Q 
with the system size reveals corrections to the topological 
value that decrease exponentially, as witnessed by the 
quality of the fits to the numerical data (symbols) that 
are shown in Fig. 7. Our fits indicate that, for all the 
three loops and for sufficiently large sizes, 



1 - e-^/« 



(11) 



where the decay rate ^ depends on the specific path. In 
particular we found « 32.35 for "Loop 1", ^2 ~ 49.98 
for "Loop 2", and ^3 w 39.10 for "Loop 3". Although a 
detailed analysis of the dependence of ^ on the chosen 
loop is computationally demanding and lies beyond the 
purposes of the present paper, we observe that ^ increases 
when the loop shrinks around the critical point. This is 
apparent by looking at Fig. 2 and noticing that, while 
"Loop 2" is entirely contained in "Loop 1", we found 
^2 > Ci- III addition, this is in agreement with the fact 
that the larger the loop, the smaller the corrections to 
the quantized value of the pumped charge due to the 
finitcness of the critical gap Ac when L is finite. 

On the other hand, when the loop docs not enclose a 
critical point, the number of pumped particles per cycle 
goes to zero in the limit i — 00, as shown in Fig. 8 for 
"Loop 4 "■ The actual asymptotic value at the thermody- 
namic limit is not exactly zero, as it is apparent from the 
data. However we point out that this discrepancy may 
be due to numerical inaccuracies (an error in Q of the 
order of £ < 0.005 is negligible on the logarithmic scale 
used for the data points in Fig. 7). 
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FIG. 8: (color online). Same as in Fig. 7, but for "Loop 4" 
(according to the notation of Fig. 2), which does not enclose 
the critical point. 



The achievement of perfect quantization (Q = or 
Q = 1) only in the thermodynamic hmit is consistent 
with the topological origin of the pumping mechanism. 
As long as L remains finite, the gap at the critical point 
Ac is never exactly zero and the topological argument 
for quantization does not hold. We have been able to 
observe quantization as long as the loop encircles the 
critical point, provided it lies at a sufficient distance from 
it. We do expect that, on choosing loops of decreasing 
amplitude, we should start observing deviations from the 
quantized value even for large (but still finite) L. For 
small loops, the number of pumped bosons per cycle, still 
retaining its geometrical origin, looses to be quantized 
unless we go to very large system sizes. Although it would 
have been interesting to show this cross-over, there are 
numerical limitations which prevent us to consider loops 
with amplitudes smaller than the "Loop 2" of Fig. 2. 
Reducing the amplitude of the loop implies a longer time 
T needed to achieve the adiabatic regime, much larger 
than T ~ 300 used to study the cycles of Fig. 2. This is 
outside our numerical capabilities. 



B. Non-adiabatic corrections 

Let us now examine in more details the non-adiabatic 
effects on the pumped charge. 

The time-resolved number of pumped particles Qo,t for 
a given finite pumping period T typically exhibits wiggles 
in time, which disappear with increasing T. In order to 
ensure the validity of the adiabatic approximation along 
the loop, one has to take a total evolution time much 
larger than the inverse square of Asm- Once this condi- 
tion is satisfied, the total number of transported bosons 
Q ~ Q{T) on a cycle is almost independent of T, thus 
promoting Q to a geometric quantity. This is apparent 
from Fig. 9, where we show Qg ( as a function of time t 
(in units of T) , for different values of the evolution period 
T and for two different sizes. The asymptotic value of the 




0.4 0.6 



FIG. 9: (color online). Time-resolved number of pumped 
particles Qo.t as a function of time, in units of the evolution 
period T. Different curves report results for different values 
of T, while the upper (lower) panel shows data for a system of 
1/ = 80 (L = 120) sites. The data have been obtained choos- 
ing "Loop 1 " and starting from the upper point in the X — V 
plane. The two right arrows indicate the final pumped charge 
Q as extracted from the large-T set of data (and eventually 
shown in Fig. 7). 
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FIG. 10: (color online). Time-resolved number of pumped 
particles Qo,t as a function of time (in units of T). We show 
results for different initial points along the path marked as 
"Loop 1", corresponding to the four extremal points on the 
V- and the A-axes. These data are for L = 80 sites and for 
an evolution time T = 200. 



total pumped charge is approached for large T, while the 
non-adiabatic wiggles tend to disappear. 

Finally we note that, once the adiabatic condition 
along the loop is satisfied, the total pumped charge does 
not depend on the initial point along the loop, as explic- 
itly shown in Fig. 10. 

We conclude this section by observing that the 
plateaus emerging in the time-resolved pumped charge 
are adiabatic, but of a non topological nature, and de- 
pend on the fine details of the instantaneous low-lying 
energy spectrum. They can change according to the spe- 
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FIG. 11: (color online). Average number of particles at each site of an open-ended chain, described by the Hamiltonian in 
Eq. (1). Different panels stand for different times during an adiabatic variation of the parameters along the path marked as 
"Loop 1". Here we chose, as initial ground state, the leftmost point of the loop corresponding to a pure MI state (see Fig. 2), 
and then follow the path in the anti-clockwise direction. The data refer to a L = 1000 chain, and to an evolution time T = 300. 
The red arrow in last panel denotes the surviving insulating region after one adiabatic loop. 



cific shape of the loop around the critical point (Fig. 9), 
as well as on the initial condition along the loop (Fig. 10). 



VII. PUMPING IN AN OPEN-ENDED CHAIN 

Quantum pumping can also be studied in an open- 
ended chain. In this case, differently from the ring geom- 
etry, the system no longer exhibits translational invari- 
ance and all the quantities defined in Sec V become site- 
dependent. In Sec. VII A we display numerical evidence 
of the fact that this configuration leads to a pumping 
mechanism which, after one loop, produces the trans- 
port of a quantized bosonic charge from one end to the 
opposite end of the chain. 

We remark that, in order to guarantee unit filling ev- 
erywhere in the chain [modelled by Eq. (1)] during the 
whole cycle, it would be necessary to couple the system to 
two (non-interacting) reservoirs attached at the two ends. 
In the absence of any direct bias between the electrodes, 
the charge would be transferred by a uniform modulation 
of the parameters A and V in the middle (interacting) 
part of the system, along the loops like those depicted in 
Fig. 2. On the other hand, in order to prevent a deple- 
tion of particles in the interacting region in favor of the 
non-interacting leads, it would be necessary to consider 
a regime of weak system-reservoir coupling, thus finding 
a trade-off between depletion and relaxation of pumped 



particles into the leads. For interacting bosonic systems, 
this regime is actually very difficult to achieve, not only in 
simulations, but even conceptually. For this reason, and 
owing to the fact that experimentally it would be much 
easier to realize an open-ended chain, we have focused on 
a system with OBC. 

A. Numerical results 

As shown in Fig. 11, one immediately recognizes that, 
during the time evolution ruled by the adiabatic variation 
of the parameters (A, V) along a closed loop (see Fig. 2), 
the average on-site density changes with time. In partic- 
ular, following "Loop 1" in the anti-clockwise direction, 
we arc able to progressively detect a net imbalance of 
bosonic particles from the right end, where we observe 
a depletion from the average unit filling, to the left end, 
where a bump in the density profile eventually appears. 
In order to probe the quantization of particle transport, 
one then has to investigate the charge Qq j pumped along 
one cycle, according to Eq. (7). 

Before discussing our results, it is important to stress 
that, for OBC, the charge transport from one edge to the 
opposite edge of the chain unavoidably drives the system 
out of equilibrium, so that, for any t > 0, the many- 
body wavefunction \'ip{t)) cannot be in the instantaneous 
ground state (ultimately this is due to the fact that, when 
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passing through the HI phase, the system faces a four- 
fold degenerate ground-state). This contrasts with the 
translationally invariant PBC setup, where the adiabatic 
approximation can be satisfied for sufficiently slow time 
variations of the Hamiltonian parameters (see Sec. IV), 
and is associated to the geometric nature of the pump- 
ing mechanism. Interestingly enough however, despite 
adiabaticity cannot be fulfilled as a global condition, it 
can be still satisfied "locally" . More precisely, this means 
that, even though the global state of the bosons might de- 
part from the corresponding instantaneous ground-state 
of the driving Hamiltonian, the reduced density matrix 
associated with a specific collection of sites remains suf- 
ficiently close to its adiabatic counterpart. Clearly, any 
deviation of the density profiles from unit filling indi- 
cates the loss of local adiabaticity. Therefore, in order to 
satisfy such condition, the instantaneous current Jj{t) 
has to be computed in the bulk of the system, where 
the particle density is still commensurate and the deple- 
tion/accumulation mechanism has not spoiled the insu- 
lating properties of the bosonic system. This requirement 
forces one to consider quite large sizes, since the region 
of deviation from the unit filling is proportional to the 
period T and involves a non negligible part of the whole 
system, in all the simulatable cases. We observed that a 
length of L = 1000 sites is adequate to our purposes'^". 

We checked the validity of the adiabatic regime locally, 
by evaluating the fidelity J- between the reduced den- 
sity matrix of two nearest neighbor sites during the time 
evolution, Pj,j+i{t) = TrL_[.jj+^{\tp{t)){il;{t)\) , and the 
corresponding one on the instantaneous ground state, 
P'j^j+iit) = Tri_{,-,+i}(|^Gs(0)(V'Gs(i)|). Omitting the 
subscripts this is given by 



J-{t) = Tr 



(12) 



The fidelity J- is a, real number which can assume any 
value in the interval [0, 1], and measures the distance be- 
tween the two quantum states. If J-" « 1, the two states 
are nearly the same, up to some global phase factor. On 
the other hand, if « 0, the two states are nearly or- 
thogonal^. In Fig. 12 we plotted J-{t) for two nearest 
neighbor sites in the central part of the chain for dif- 
ferent values of the period T (top panel), and for the 
largest period as a function of the position j of the link 
in the chain (bottom panel). For a total evolution time 
T = 300, T{t) w 1 along the whole loop, thus provid- 
ing evidence of adiabaticity, in the central region (top 
panel of Fig. 12). By departing from the middle, even 
at T = 300 adiabaticity is lost (see the bottom panel of 
Fig. 12). The results for the pumped charge reported 
below are taken in the region close to j = 500, where 
adiabaticity is ensured during the whole cycle. 

Following this prescription, we averaged the current 
over the sites marked by the red arrow in the last panel 
of Fig. 11, where we saw deviations from n = 1 of up to 
10~^, and the local adiabatic condition is satisfied. The 
results are shown in Fig. 13, for three different values 
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FIG. 12: (color online). Two-site nearest-neighbor fidelity 
between the instantaneous state during the time evolution 
and the corresponding adiabatic ground state for a chain with 
L = 1000. (Top panel) The fidelity for two neighbouring sites 
in the centre of the chain is analysed as a function of the 
period T. For T — 300 the deviation from adiabaticity is 
below 10~^ during the whole evolution. (Bottom panel) The 
fidelity is analysed as a function of the position of the link 
in the chain for the largest period (T = 300). The observed 
deviations close to the end of the chains are in correspondence 
with the deviations from unit filling visible in Fig. 11. 



of the evolution period T along "Loop 1" in the anti- 
clockwise direction. For the maximum T wc considered 
(T ~ 300), we observe a remarkably high level of ac- 
curacy in the quantization of the pumped charge. For 
lower values of T, perfect quantization is lost and wiggles 
signaling adiabaticity losses start to emerge, in analogy 
with the ring geometry already discussed in Sec. VI B. 
To check the consistence of this result, we also computed 
the total number of particles on the left and on the right 
of the insulating region, finding that deviations from half 
filling correspond, at the end of the cycle, to having re- 
spectively one particle more and one particle less. 

Finally we simulated the pumping protocol by invert- 
ing the travelling direction along the loop, and found 
a scenario qualitatively similar to that of Fig. 11, with 
the depleted and the filled edges inverted. The corre- 
sponding pumped charge behaves analogously, despite 
the pumping mechanism is working in the opposite direc- 
tion. We point out that the two directions are not sym- 



10 




— T = 


100 


— T = 


200 


— T = 


300 



0.4 0.6 

t/T 

FIG. 13: (color online). Time-resolved number of bosons 
which are pumped from the right to the left side of the open- 
ended chain, as obtained from Eq. (7). The various continuous 
lines refer to a system of L = 1000 sites, for three different 
evolution periods T along "Loop 1", starting from the pure 
MI point and going in the anti-clockwise direction. The in- 
set displays the instantaneous current (t) averaged over the 
insulating region highlighted in Fig. 11, for T = 300. 



metric, while only the total number of pumped bosons is 
quantized in the two cases. 



rameter space leads to a quantized number of pumped 
particles with topological protection. 

A possible realization has been proposed by Berg et 
al}^ and consists of an extended Bose- Hubbard model in 
one dimension in the presence of an inversion-symmetry- 
breaking term. Here we have investigated this model 
by means of time-dependent density matrix renormal- 
ization group simulations, giving numerical support to 
their original prediction. We have computed the pump- 
ing properties of the system as a function of the system 
size and shown that the scaling is consistent with a num- 
ber of pumped bosons per cycle equal to one or zero, 
depending on whether the loop performed in parameter 
space encircles or not the critical point. In addition, we 
have provided evidence that the pumping mechanism is 
indeed adiabatic and that it displays topological protec- 
tion by showing that, in the thermodynamic limit, the 
number of pumped particles is independent of both the 
total evolution time T and of the specific path performed 
in parameter space. We considered two different config- 
urations: a ring geometry, such that a persistent current 
arises as a consequence of the pumping mechanism, and 
an open chain, where quantum pumping is manifested 
in a depletion of charge at one end of the chain with a 
corresponding accumulation at the opposite end. 



VIII. CONCLUSIONS 
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